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Abstract 

One of the largest theoretical uncertainties assigned to the strong coupling con- 
stant a s as determined from hadronic tau decays stems from the differences in the 
results for Fixed Order Perturbation Theory (FOPT), Contour Improved Pertur- 
bation Theory (CIPT) and Renormalon Chain Perturbation Theory (RCPT). It is 
often argued that the three methods differ in the treatment of higher orders only and 
therefore the full difference should be treated as theoretical error. Recently other 
arguments either in favor of FOPT, CIPT or RCPT have been given, but none of 
those is able to combine all three to a single value in the strong coupling constant. 
In this note I will show that FOPT alone has a much larger uncertainty than pre- 
viously assumed and therefore agrees within error with CIPT. Furthermore a more 
appropriate matching of the different schemes used in RCPT reduces the difference 
to the CIPT result by a factor of 6. Together with recently published results for 
the 4th order term K4 this reduces the theoretical error on a s by a factor of 2.5 
compared to the previously assumed spread of the three perturbative approaches. 

1 Introduction 

Hadronic decays of the r lepton are among the most actively studied fields in QCD. 
The unique situation of a small mass scale and still small non-perturbative contributions 
allow for a very precise determination of the strong coupling constant a s [1-4] (For recent 
reviews see [5-10]). The ratio of the hadronic decay width of the r and its leptonic width 
can be written as 

R T = 35*EW (|Kid| 2 + |Kis| 2 ) (1 + #EW + ^P crt + ^non-port) , (1) 

where S E w = 1.0198±0.0006 [11, 12] and 5' EW = 0.0010±0.0010 [13] are small electroweak 
corrections, <5 non _ pert denotes a 0(few%) non-perturbative correction and <5 pert is the per- 
turbative prediction. Neglecting the masses of the quarks (as is a good approximation for 



1 



the non-strange decay width of the r) the perturbative part is given by 



\s\=mi 

with the known coefficients [14-19] 

K = K 1 = l, 
299 

K2 = |f-9C(3), 

kF = 5805T_rr9 c(3) + ?c(5); (3) 

Ms 78631453 4185 2 1704247 34165 1995 

^ 4 = ^0736^ + ~ C(3) " ^32~ C(3) + ^6~ C(5) " "TF C(7) ' 

where K 3 and K4 are scheme dependent and given here in the MS-scheme and £(n) denotes 
the Riemann zeta function. The 4-th order term K^- s ~ 49.0757 deviates substantially 
from previous estimates and partial calculations of that coefficient 7^P artial — 27 ± 16 [20, 
21]. The fifth-order term has been estimated to K 5 ~ 275 in [19], but the large deviation 
of the exact K 4 from it's prediction suggests that a 100% error on K 5 is realistic. For 
the purpose of evaluating differences stemming from the 5th and higher orders I'll use 
K 5 = 400 ± 400 in this note. 

The methods FOPT and CIPT [4] differ in the way (2) is calculated. In the CIPT 
approach the /^-function is used to get numerical solutions for a s (— s) in the complex 
s-plane by starting with a s {m 2 ). The integrand is thus calculated in small steps on the 
circle \s\ — m 2 and the sum of all pieces gives the total integral. 

For the FOPT method the /3-function and its derivatives are Taylor expanded in s 
around sq = m^. which leads to a power series representation of a s (—s) in powers of 
a s (rri^). The series is truncated at the desired order (here the 5th) in the strong coupling 
and inserted in the integral which becomes solveable now. The usual FOPT result reads: 

<W = ^ + 5.2023^1 + 26.366^^ + (4) 

127.08^^ + (K 5 + 307.78)^^ + 0(a s 6 ). 

As is demonstrated in [19,21] the fourth and fifth order terms contribute very little to 
the perturbative part and the difference between the FOPT and the CIPT result is much 
larger than the contributions from these terms even if generous errors are used for K 5 . 
Taking a s {m*) = 0.35 and K 5 = 400 as reference values we could first calculate <5 pcrt from 
the CIPT approach and extract a s {m^) again using FOPT: 

6 T (o,(m?) = 0.35) = 0.21179, (5) 

a s (m 2 T ) FOPT \5 pcrt = 0.21179) = 0.32543. (6) 

The deviation of either value from their mean is with Aa s = ±0.012 almost twice as 
large as the uncertainty due to higher orders Aa sAK5 = 0.007. The reason for this large 
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difference is the choice of the point on the circle \s\ — in the complex s-plane around 
which the /3-function and its derivatives are Taylor expanded to approximate the strong 
coupling on the circle. In the following section the FOPT formalism will be generalized 
to allow for other choices. 



2 Generalized FOPT 

The starting point is the perturbative expansion of the (3 function, which is given by 

P(a s ) = ^ = -p al(s) - ftofr) - (3 2 a%s) - (3 3 a 5 s (s) - . . . , (7) 

where a B (s) = cc s (s)/(47r). The first two terms in the /3-function [22-26] for m quark 
flavors, 

2 

A> = 11 - g^f^ 
38 

Pi = 102 -—n f , 

are universal at leading twist whereas the higher order terms are scheme dependent. In 
the MS scheme the first two scheme dependent coefficients are known [27-30]: 

„t™ 2857 5033 325 2 

ft = — "TT n ' + lA n " 

„ m 149753 w , /1078361 6508 ^ A 

/50065 6472 w A , 1093 , 

The Taylor expansion of the evolution equation (7) around s reads up to the fifth 
order in a s : 

gg(g) = u s (so) l^ ln s ( a s {s Q ) V | 

7T 7T 4 So \ 7T / 

16 V s s J \ n J 

128 V S $0 S J \ IX J 

-i-(6/3 4 ln 4 A_26/3 2 Aln 3 ^+ 
1536 V s s 

^) +OK(s ) 6 ). (9) 

It should be noted that eq. (9) is strictly speaking not a Taylor approximation since the 
truncation occurs at a certain power of a s and not at a certain power in the expansion 
variable In (s/s ). 
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Fig. 1: Quality of the Taylor expansion ofa s {m^ exp(^)). The two plots show the absolute 
value \oi s \ on the complex circle s = m^exp(zyj). The left plot shows with long- dashed, 
dash-dotted, narrow dotted and wide dotted lines the Taylor (FOPT) expansion up to 5th, 
4th, 3rd, and 2nd order, respectively. The J^-loop result for the numerically solved (3- 
function ( CIPT) is drawn as solid line for comparison. The right plot shows with dahsed, 
dotted, and dash-dotted lines the numerical solutions for 3, 2 and 1 loop ^-functions, 
respectively. Again the same 4-loop result as in the left plot is shown as a solid line. The 
reference value a s (m^) = 0.35 was used for all curves. 



Since both the integrand in eq. (2) and eq. (9) are power series in a s it is interesting 
to compare the magnitudes of the coefficients in these series. The largest values in eq. (9) 
are obtained at s = Soexp(— in) where the |c n | (the magnitude of the coefficient in front of 
(a s /7r) n ) read 1, 7.07, 51.52, 390.8, 3023.85. At s = s exp(— in/2) (the average distance 
from So) the | 1, 3.53, 13.98, 62.33, 275.16. These numbers grow much (slightly) 

faster at — n (— 7r/2) than the K n , where (starting with K\) we have 1, 1.64, 6.37, 49.08, 
~ 275. Therefore it is conceivable that the nature of the Taylor expansion of a s dominates 
the uncertainty of the FOPT result and not missing higher order K n terms. To illustrate 
this the Taylor expansion of a s is modified in the following. 

Figure 1 shows that the deviation of the Taylor expanded ct s from the numerically 
solved a s grows with the distance of s = soexp(^) from the chosen development point 
So = vr? T . The CIPT results for 1 to 4-loop treatment on the circle are also shown. Starting 
at 2-loop level the CIPT results are almost indistinguishable while the FOPT deviations 
remain large even at 5th order! Thus the correct treatment of the logarithms and not 
higher orders in a s are the dominant source of uncertainty. 

Furthermore the choice so = vri^ for FOPT is completely arbitrary. This becomes even 
more obvious in the light of the usual procedure of comparing different a s measurements 
by evolving them to the Z°-mass with a numerically solved /3-function. If we were to evolve 
the a s from tau-decays to the Z°-mass in a FOPT-like manner with equation (9) in just 3 
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steps (with the quark flavor transitions to rif — 4 and rif — 5 at m T and mb, respectively) 
the numerical value of a s (m 2 z ) (for a s (m 2 T ) = 0.35) would be 0.042, 0.206, 0.064, 0.160, 
for the 1, 2, 3, 4 loop beta function, respectively. Compared to the usual procedure of 
contour- improved evolution of the coupling (i.e. using eq.(9) in small steps) which yields 
0.128, 0.123, 0.122, 0.122, for 1, 2, 3, 4 loop, respectively, the 3 step FOPT solutions 
converge very slowly, giving alternatingly lower and higher estimates of the coupling as 
more and more orders in a s are considered and still at 5th order resulting in a numerical 
value that is far below the exact result. The FOPT terms up to the 5th order do not 
compensate for the neglected large logarithms. 

It is therefore natural to generalize this in the case of the r and first evolve a s (m 2 ) 
to a s (m 2 exp(i<£> )) with the numerically solved /3-function and derive the Taylor series of 
<5 pert around this new point. The integral (2) can in fact be split in two pieces around 
ifi and — ifi since the strong coupling at — y?o is just the complex conjugate of the strong 
coupling at y9 : 

a s (m 2 exp(-i</? )) = a s (m 2 T exp(iy9 ))* (10) 
The resulting S pert up to the fifth order reads: 
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or numerically: 

<W = a + 0.8488 6 + (-4.5 y? + 0.8082) a b + 

(1.9099 p + 5.2023) (a 2 - 6 2 ) + 

(-5.0625 ip 2 + 5.2138 y? + 26.366) (a 3 -3a b 2 ) + 

(4.2972 v?o + 27.410 y9 + 12.356) (b 3 - 3 a 2 6) + (12) 
(-9.6687 (pi - 101.51 v?o - 71.629 v?o + 127.08) (a 4 - 6 a 2 6 2 + 6 4 ) + 
(45.563 ifl - 100.94 y? 2 - 918.59 y9 - 521.11) (a 3 b - a b 3 ) + 
(25.629 v?o - 92.897 v?o - 1220.5 y? 2 , - 1272.5 y? + 

X 5 + 307.78) (a 5 - 10 a 3 b 2 + 5 a 6 4 ) + 
(21.755 v?o + 324.78 y? 3 , + 271.83 y? 2 , - 1612.0 y? + 

0.8488 K 5 - 1413.5) (6 5 + 5 a 4 6 - 10 a 2 b 3 ), 

with </?o € [— 7r, 0], and a s (m 2 exp(iy2 ))/7r = a + ib. Three points should be noted about 
equation (11): 

1. it resembles the usual FOPT result for </?o = and 6 = 0; 

2. <5 P ert remains real for all choices of a s and ipo] 

3. inserting the Taylor expanded a s (y?o) in eq. (11) and Taylor expanding the resulting 
<5 P ert again around a s ((p = 0) leads also to the usual FOPT result. 

The last point demonstrates that FOPT can be generalized only if the 'exact' value for 
a s ((po) is used in the expansion. 

Figure 2 shows 5 pcit as a function of y? with a s (m 2 ) = 0.35 and K 5 = 400,0,800 as 
reference values. The consequences of the generalized FOPT solution are discussed in the 
following section. 



3 Discussion of the generalized FOPT solution 

As can be seen from figure 2 the FOPT result depends largely on the choice of (f . The 
FOPT curves intersect with the CIPT curves around ipo — —1 but span over a much larger 
range of <5 pe rt values. Compared to the uncertainty from the neglected higher orders this 
intrinsic error is 4 times larger as none of the choices for ip should be excluded. The 
default choice of <p — leads to the largest possible value of <5 pcrt and therefore a s from 
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Fig. 2: t5 P erf as function of the development point (po. The solid, medium-dashed and short- 
dashed lines show the CIPT result to 5th order for K 5 = 400, 0, and 800, respectively. 
The long-dashed, dash-dotted and dotted lines show the generalized FOPT result to 5th 
order for K5 = 400, 0, and 800, respectively. The reference value a s (m^) = 0.35 was used 
for all curves. 



FOPT used to be smaller than from CIPT. The deviation can however not be attributed 
to higher order terms in the series of d~ peT t- Instead the extraction of a s with FOPT 
should use the average of the two extremes a s ((Po = 0) and a s (<^o = ~ and half of 
their difference as additional theoretical error. Consequently the most accurate way for 
the determination of a s from r decays is the CIPT approach. There is no reason to add 
the same error to the CIPT result as it does not depend on the choice of ipQ. Also since 
the FOPT result agrees within its own error with CIPT there is no discrepancy anymore 
between results with these two approaches. In fact the CIPT result is what FOPT would 
converge to for n — > oo if n equidistant points on the circle s = Soexp(i^) would be used 
in the expansion. The case n = 1 could therefore be regarded as an approximation for 
CIPT and the choice <po = is just one of the many possible choices for n — 1. 
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4 Renormalon Chains 



The third theory often used in evaluating a s from r decays uses so called 'Renormalon 
Chains' [31-34] and re-sums the f3o parts of S pert to all orders in a s . In [32] the result 
obtained from this re-summation is corrected by the known FOPT terms by first sub- 
tracting the large-/5 part up to the desired order of FOPT and then adding the FOPT 
part. Thus 5 pcTt for the Renormalon Chain Perturbation Theory (RCPT) can be written 
as 

rRCPT _ r _ rFOPT , rFOPT /-, o\ 

"pert — "renormalon "large-/3 ' "pert i \ 10 ) 

where the three terms in the sum refer to the renormalon chain result, the large-/3o re- 
summed result up to the order used in FOPT, and the FOPT result, respectively. As is 
pointed out in [32] the renormalon chain re-summation includes parts of the CIPT re- 
summation, namely the terms (— /5 /4 ln(s/s )) n which are part of the coefficient in front 
of (a s /7r) n in eq. (9). Therefore the CIPT result can not be used instead of the FOPT 
result in eq.(13). Still, for the FOPT correction and the fixed order large-/3o correction 
the same arbitrariness of the choice of ip as discussed in the first part of this note exists, 
as long as both the FOPT term and the fixed order large-/9 term are expanded around 
the same ip . Therefore the variation of 5* — £ia? g c-/3 w ^h V^o is a source of uncertainty 
in the RCPT approach. The generalized o\a?ge-/3 can be derived from eq. (11) by setting 
f3 n = for n > and replacing the K n with [3$™ ^K n , which are given up to n = 4 in [32] 
and up to n = 12 in [10]. Numerically o^^e-A U P to the fifth order in a s is given by: 

5™ PT = a + 0.8488 b + (-4.5 ip + 0.6668) a b + 

(1.9099 p + 5.1190) (a 2 - b 2 ) + 

(-5.0625 y? 2 + 1.5002 <p + 28.779) (a 3 -3a b 2 ) + 

(4.2972 ip 2 + 23.035 v? + 2.5067) (6 3 - 3 a 2 b) + (14) 
(-9.6687 v?q - 77.745 ip 2 - 16.920 v?o + 156.67) (a 4 - 6 a 2 b 2 + b 4 ) + 
(45.563 ip 3 - 20.253 y? 2 - 777.03 (p - 433.69) (a 3 b - a b 3 ) + 
(25.629 v? 4 , - 15.189 v?o - 874.16 <p 2 - 975.80 v? + 

900.78) (a 5 - 10a 3 b 2 + 5a6 4 ) + 
(21.755 v?o + 233.23 y? 3 + 76.141 y? 2 - 1410.1 y - 

6 1 5.93) (6 5 + 5 a 4 6- 10 a 2 6 3 ), 

with ip G [— 7r, 0], and a s (m 2 exp(i<p )) / V = a + ib as in eq. (12). The difference 5p°f T — 
<5iarge-/? is n °t as sensitive to the choice of ipo as 5p°f T alone and roughly halfs the 
associated uncertainty in a s . 

Figure 3 shows the RCPT result using the (^renormalon as in [32] but the modified 5* °f T — 
<5iarge-/? from eqs. (12,14) to correct the result up to the fifth order in a s . The reference 
value of a s {m 2 ) = 0.35 is used again for all curves. It is clear from the figure that RCPT 
would still require a much smaller a s ~ 0.31 compared to CIPT even with the modified 
corrections in the fixed order parts. 

This observation relies however on the fact that a s in the renormalon part and the 
fixed order part of eq. (13) refers to the same quantity. This is probably not the case. The 
renormalon part in [32] is derived from the one-loop coupling in the so-called V scheme, 
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Fig. 3: 5 per t as function of the development point tpo. The solid, medium-dashed and short- 
dashed lines show the CIPT result to 5th order for K 5 = 400, 0, and 800, respectively. 
The long-dashed, dash-dotted and dotted lines show the RCPT result with 0-loop matching 
and generalized FOPT correction to 5th order for = 400, 0, and 800, respectively. The 
reference value a s {m 2 ) = 0.35 was used for all curves. 



a^(fi 2 ) which is matched on the 0-loop level to a^ s (exp(— 5/3)/i 2 ). The problem therefore 
is that we have a coupling constant on the 3-loop level 1 in the FOPT parts, but treat it as 
a one-loop coupling in the renormalon parts. A possible solution would be to use 2-loop 
matching to go from the MS-scheme to the V scheme which is given by [35, 36] 



(//) _a 8 MS (e- 5 V) _ ^ MS (e°- 422 V 2 )V 7 , 72£16 ( a ™ (e a42 V ) , , (15) 



7T 7T \ 7T / \ 7T 

Figure 4 shows again the RCPT result as before but with the 2-loop matching for a 



The reference value of a s (m 2 ) = 0.35 is used for all RCPT and CIPT curves. The large 



1 Since enters only in the 5th order in R T a s is effectively used as a 3-loop coupling constant in 
FOPT which goes up to the 4th (plus estimated 5th) order in a s 
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Fig. 4: 5p er t as function of the development point (fo. The solid, medium- dashed and short- 
dashed lines show the CIPT result to 5th order for K 5 = 400, 0, and 800, respectively. 
The long-dashed, dash-dotted and dotted lines show the RCPT result with 2-loop matching 
and generalized FOPT correction to 5th order for K§ = 400, 0, and 800, respectively. The 
reference value a s (m^) = 0.35 was used for all curves. 



overlap of the CIPT and RCPT curves shows that the differences in the deduced strong 
couplings from both theories are much smaller than previously assumed. 



5 Discussion of the modified RCPT solution 

Unlike in the case of the generalized FOPT the RCPT solution can not be regarded as the 
first iteration of a contour improved result as this would result in inconsistent definitions 
of a s on the circle s = \m%\. Therefore the magnitude of the spread of a s values obtained 
from a fixed 5 pert and a fixed K 5 could be regarded as induced by higher order terms. 
This spread has furthermore the same magnitude as the difference between CIPT and 
RCPT at <£> = 0, which at the same time shrunk by a factor of six by using the 2-loop 
matching instead of 0-loop matching for the transition from the V scheme to the MS 



11 



0.38 

B 

8 

0.36 



0.32 



0.3 







a ciPT k _ o 

s 5 — 

af PT K- = 400 
af PT K. = 800 




* 

* 


\ 

S 


— •■ 








* 


■ ■■■■■ IliLI ■ ■ 













St! 


•£•£! 




si****" 










a RCPT K 












•••••• 


a« CPT K 
a« CPT K 


5= 

5 = 400 
5 = 800 









■2.5 



-1.5 



-0.5 





% 



Fig. 5: a s (m%) as function of the development point (fo. The solid, medium- dashed and 
short-dashed lines show the CIPT result to 5th order for = 400, 0, and 800, respec- 
tively. The long-dashed, dash-dotted and dotted lines show the RCPT result with 2-loop 
matching and generalized FOPT correction to 5th order for K$ = 400, 0, and 800, respec- 
tively. All curves are obtained with 8 pert = 0.21179 as reference value. 



scheme. Averaging over all <^ = values for RCPT leads to the same numerical value 
for a s as CIPT at the central choice for K§, while the RCPT results for Ks = 0, 800 stay 
much closer to the central a s compared to the corresponding CIPT values, showing that 
the large (3q re-summation reduces the impact of higher order terms. 

Recently the situation about the influence of higher-order corrections has been revis- 
ited in [10], where the authors study the renormalon structure in different models and fix 
the Borel transform of the Adler function such that the known terms up to 4th order are 
reproduced. Using this matched Adler function for R T and comparing the full re-summed 
result with standard FOPT shows again good agreement, while the distance to CIPT is 
large. But as discussed in section 2, the generalized FOPT solution would on average 
reproduce the CIPT result and therefore the model of [10] would also deviate from the 
approach advocated in section 2. 
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The power corrections to R T have been re-examined in [37] where duality violation 
parts [38,39] of the order of 0.01 (but found to be negligible in [9]) and tachyonic mass 
corrections from the gluon are considered. It is argued in [37] that the difference between 
the Borel-sum and the truncated series at 4th order in the large-/9o limit can be regarded 
as non-standard dimension 2 power corrections to R T . Since these corrections are of the 
order 0.04 for CIPT and 0.02 for FOPT, respectively, they dominate over the duality 
violation effects and if taken at face value eliminate the difference between FOPT and 
CIPT. However, it should be noted that the Borel transform D(b) of the Adler function 
does not have a pole at b = 1 and therefore no renormalon ambiguity of dimension 2, 
which makes it difficult to associate a dimension 2 correction to the observed difference. 

Here I'll concentrate on the perturbative parts only and neglect any non-standard 
power correction. Figures 2 and 4 show that there is no large cancellation mechanism 
which would prevent the FOPT or the RCPT result from depending strongly on the 
arbitrary choice of the development point ipo. Including this arbitrary choice in the 
uncertainty estimate shows that CIPT provides the most accurate estimate. It is however 
re-assuring that correcting the large-/?o re-summed result with the known fixed order terms 
up to 5th order reduces this dependency by 50% leaving the remaining 50% to the exact 
logarithms for (5 n with n > 0. 

Figure 5 shows a numerical example for S pert = 0.21179 comparing the values of a s for 
K 5 = 400, 0, 800 from fits to CIPT and RCPT as a function of tp Q . The CIPT numbers 
are: 

a s (m 2 T , S pert = 0.21179, K 5 = 0) CIPT = 0.3572, 

a s (m 2 T ,S pert = 0.21179,^ = 400) CIPT = 0.35, (16) 

a s (m 2 T , S pcrt = 0.21179, K 5 = 800) CIPT = 0.3432, 

and for RCPT the result is: 

a s (m 2 T , 5 pcrt = 0.21179, K 5 = 0) RCPT = 0.3519 ± 0.0097, 

a s {m 2 T , 5 pcrt = 0.21179, K 5 = 400) RCPT = 0.3498 ± 0.0063, (17) 

a s (m 2 T , 5 pcTt = 0.21179, K b = 800) RCPT = 0.3480 ± 0.0034, 

where the errors are given by the RMS of the a s values over the range — n < ip < 0. Now 
RCPT and CIPT agree on the central value of a s and both give similar estimates for the 
uncertainty due to (different) neglected higher order terms: 

a s «, 5 pcrt = 0.21179) = 0.3499 ± 0.00721°;°°^. (18) 

where the first error is due to K 5 ~ 400 ± 400 and the second due to the variation of the 
renormalization scale 0.4 < /i 2 /m 2 < 1.6. 

6 Numerical analysis 

Using the same numerical value for S peTt = 0.2042 ± 0.0038 cxp ± 0.0033 non _ per t as obtained 
in [9] and used in [10], where the first error is the experimental one, dominated by the 
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non-strange hadronic decay ratio of the r, R t ,v+a and the second is due to the non- 
perturbative and quark-mass corrections, the results for CIPT and generalized FOPT 
and RCPT read: 

af PT (m?) = 0.3406 ± 0.0047 exp ± 0.0041 non _ pert ± 0.0066^, 

«' OP Vr) = 0.3535 ±0.0061 exp ±0.0053 non _ pert ±0.0208^lHooi^ 6 , (19) 

af CPT (m 2 T ) = 0.3440 ± 0.0030 cxp ± 0.0026 non _ pcrt ± 0.0061 vo ± 0.0019^ 5 ., 

All three results agree within the error due to tpo which is very large for FOPT but 
moderate in case of RCPT. The difference between CIPT and RCPT is of the same size 
as the error due to </? for RCPT and the average between both values (and conservatively 
assigning the larger of the two results errors to the average) leads to: 

a s (m 2 T ) = 0.3423 ± 0.005 cxp ± 0.007 Ai , 5 ± 0.004 non _ pcrt i°;°°^, ( 2 0) 

where the fourth error is due to the variation of the renormalization scale. The total 
theoretical error (including the non-perturbative part) is with ±0.008 only marginally 
larger than the experimental error. 

Evolving a s given by eq. (20) from m r = 1.7768 GcV to m z o = 91.1876 with m c {m c ) = 
1.27lo:?iGeV and m b (m b ) = 4.20±o;o?GeV [40] with the flavor thresholds to n f = 4 and 
71/ = 5 at m T and mb, respectively, gives: 

a s {m% ) = 0.1213 ± 0.0006 cxp ± 0.0008 A ^ 5 ± 0.0004™ n _ pcrt t°;°°°^ ± 0.0002 cv , (21) 

where the last error is the evolution uncertainty due to the variation of the thresholds 
m q < m t hrcsh < 2m g and the quark masses itself within their respective errors. 



7 Conclusions 

Modifying the usual Taylor expansion in fixed order perturbation theory by allowing 
starting points other than a s (m 2 ) on the complex circle \s\ = m 2 . reveals a larger intrinsic 
uncertainty of FOPT than previously assumed. Giving equal weight to all possible choices 
and averaging over the different values of a s FOPT so obtained brings the FOPT result in 
agreement with CIPT. Since CIPT does not bear this additional intrinsic uncertainty the 
CIPT solution should be preferred. The large-/3 re-summed result can be modified in its 
fixed order parts with a similar approach. Here the variation due to the starting point 
of the strong coupling on the complex circle alone does no account for the difference to 
CIPT. But applying 2-loop matching instead of 0-loop matching to combine the large-/3o 
parts with the known fixed order parts up to cancels the difference of RCPT and CIPT. 
In this way all three perturbative approaches finally agree on the central value of a s from 
the r. The final result from eq. (21) is 

a s (m 2 z0 ) = 0.1213 ± 0.0006 cxp ± 0.0010 th eo, (22) 

where the first error is experimental and the second theoretical including power correc- 
tions. This confirms with different theoretical arguments the large a s obtained in [9] and 
is not compatible with the numerically lower values in [10,37,41,42]. 
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